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We propose a new mechanism for explaining the observed asymmetry between 
matter and antimatter, based on nonpcrturbative physics at the QCD scale. Our 
mechanism is a charge separation scenario, making use of domain walls separat- 
ing the recently discovered long-lived metastable vacua from the lowest energy 
vacuum. The walls acquire a fractional negative baryon charge, leaving behind 
a compensating positive baryon charge in the bulk. The regions of metastable 
vacuum bounded by walls ("B-shells") will contribute to the dark matter of the 
Universe. 



1 Introduction 

The origin of the asymmetry between baryons and antibaryons, and, more 
specifically the origin of the observed baryon to entropy ratio tib/s ~ 10~ 10 
(tib being the net baryon number density and s the entropy density) remains 
a mystery and one of the main challenges for particle-cosmology. In order 
to explain this number from symmetric initial conditions in the very early 
.Universe, it is generally assumed that three criteria, first laid down by Sakharov 
El must be satified: 

• Baryon number violating processes exist. 

• These processes involve C and CP violation. 

• The processes take place out of thermal equilibrium. 

Here we report on a recent proposal! that baryogenesis may be realized at 
the QCD phase transition. This scenario is based on the existence of domain 
walls separating the recently discovered u long-lived metastable vacua of low- 
energy QCD from the stable vacuum. The walls acquire a negative fractional 
baryon charge, leaving behind a compensating positive baryon charge in the 
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bulk. In this sense, our proposal is a charge separation rather than a charge 
generation mechanism. 

It is well known that topological defects can play an important role in 
baryogenesis. This was already realized □ in the context of GUT scale baryoge- 
nesis u. In grand unified theoreis (GUT), baryogenesis can occur at or immedi- 
ately after the symmetry breaking phase transition via the out-of-equilibrium 
decay of the superheavy Higgs and gauge particles X and a perturbative 
process. Defects provide an alternative mechanism. If during the symmetry 
breaking phase transition defects are formed, then a substantial fraction of the 
energy is trapped in these defects in the form of topological field configura- 
tions of X and A^. Upon the decay of the defects, the energy is released as X 
and A^ quanta which subsequently decay, producing a net baryon asymmetry. 
Note that defects below the phase transition represent out-of-equilibrium field 
configurations, thus ensuring that the third Sakharov criterium is satisfied. 

As was realized by Kuzmin et al. i, any net baryon asymmetry produced 
at very high energies can be erased by n baryon number violating nonpertur- 
bative processes (sphaleron transitions 0) which are unsuppressed above the 
electroweak scale. Hence, a lot of attention turned to electroweak baryogen- 
esis, the attempt to re-generate a nonvanishing ub/s by means of sphaleron 
processes below the electroweak scale, when they are out of equilibrium (see 
e.g. lj for recent reviews). 

Topological defects may also play a role in electroweak baryogenesis i. If 
new physics just above the electroweak scale generates topological defects in 
the cores of which the electroweak symmetry is unbroken, then these defects 
can mediate baryogenesis below the electroweak scale. The defects are out-of- 
equilibrium field configurations. In their cores, the sphaleron transitions may 
be unsuppressed, and, typically, C and CP violation is enhanced in the defect 
walls, thus demonstrating that all of the Sakharov criteria are satisfied. 

Detailed studies (see e.g. 0) , however, have shown that without introduc- 
ing new physics (e.g. supersymmetry with Higgs and stop masses carefully 
chosen to lie within narrow intervals), electroweak baryogenesis is too weak 
to be able to generate the observed value of ris/s ~ 10~ 9 . This criticism 
applies in particular to string-mediated electroweak baryogenesistU. Thus, at 
the present time the origin of the observed baryon to entropy ratio remains a 
mystery. 

It is therefore of interest to explore the possibility that the baryon asymme- 
try may have been generated at the QCD scale via nonperturbative processes, 
without the need to introduce any new physics beyond the standard model, 
except for a solution of the strong CP problem. Since baryon number is glob- 
ally conserved in QCD, the only way to produce a baryon asymmetry is via 
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charge separation. It is to a discussion of such a mechanism to which we now 
turn. 



2 QCD Domain Walls 

Crucial for our scenario is the existence of QCD domain walls, a consequence 
of the recent improved understanding of the vacuum structure of QCD. 

Let us first consider a SU(N) gauge theory in the absence of fermions. 
Since the gauge configuration space has nontrivial topology, there are discrete 
states |n > which minimize the energy in each of the subspaces of field con- 
figuration with winding number n (an integer), and from which in turn the 9 
vacua \9 > can be constructed: 

\9>=J2e m9 \n> . (1) 

By construction, physical quantities, in particular the ground state energy, 
must be 27r-periodic in 9. The effects of 9 ^ can be recast into an additional 
term in the QCD Lagrangian: 

Cqcd = ^Tr V" + JL^TrF^Fap , (2) 

where F^ v is the gauge field strength tensor and g is the coupling constant. 

In the large N limit, the combination A = g 2 N must be held constant. 
Hence, in this limit the vacuum energy as a function of N and 9 must scale as 

E(0,N) = N 2 h{9/N) (3) 



where h is a continuous function. It is also knownE3lL3 that E(9, N) is nontrivial 
(and proportional to 9 2 ) for small values of 9. It is difficult to reconcile this 
with (||) and with the 27r-periodicity of E(9) unless E(9) has a multi-branch 
structure 

E{9) = N 2 min k h((9 + 2irk)/N), (4) 

where h is a smooth function. Such a multi-branch structure was first pro- 
posed for supersymmetric QCD F 4 H 15 I. In a functional integral approach, the 
prescription corresponds to summation over all branches in a multi-valued ef- 
fective Lagrangian. In the thermodynamic limit, only the branch with the 
lowest energy contributes. 

Let us now include fermions. In the low energy limit, only the pions ir a 
and the rf meson contribute. They can be described in terms of the matrix 

U = exp[iV2—^ +i-=f , (5) 
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where A a are the Gell-Mann matrices and f„ (jy) is the pion (77') coupling 
constant. 

According to the anomalous Ward identities, for massless quarks (m q = 0) 
the ground state energy as a function of 9 and U can only depend on the 
combination 9 — iTrlogU 

E{9, U) = E{9 - iTrlogU) . (6) 

From (^) we can immediately derive the form of the effective potential for U 
for a fixed value of 9: 

V eff {U) = E(9,U) = E {9 - iTrlogU) , (7) 

which inherits the multi-branch structure of the function E{9) of the pure 
gauge theory. In particular, there are distinct ground states. In the chiral 
limit, their energies are degenerate, but for finite quark masses the degeneracy 
is softly broken. The energy density barrier between neighboring vacua is of the 
order Aq CD which is much larger than the energy density difference between 
the minima which is0 dp ~ m q h^ CD (Aqcd is the QCD symmetry breaking 
scale, about 200MeV). The walls are described by a tension a which is of the 
order A 3 QCD . 

The presence of degenerate minima of V e ff(U) leads to the existence of 
domain walls separating regions in space where U has relaxed into different 
minima. Such domain walls will inevitably be formed O during the QCD 
phase transition. 

As first realized inH, the presence of QCD domain walls implies that the 
second and third Sakharov criteria are automatically satisfied. Note that per- 
turbative QCD processes in the different minima are the same modulo the value 
of the strong CP parameter which is shifted by A9j — iTrlogU j in the j-th 
minimum Uj . In order to avoid the strong CP problem, the effective value 9 e f f 
must be close to zero in the global minimum at the present time. In this case, 
8 e ff will be of the order 1 in the meta-stable vacua. Hence, there is (almost) 
maximal CP violation across the domain walls. Note that no new physics is 
required in order to generate a large amount of CP violation. As stressed in the 
Introduction, the domain walls are out-of-equilibrium configurations. Hence, 
Sakharov's second and third criteria are satisfied. 

3 Induced Charge on the Domain Wall 

It has been known for a long time 111 that solitons can acquire fractional 
fcrmionic charges. The prototypical example is a 1 + 1-dimensional theory 
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with fcrmions coupling to a real scalar field with a double well potential via a 
Yukawa coupling term. In the background of the kink solution for the scalar 
field, the effective Lagrangian for the fermions ifi is 

£ 2 = ^{idjT 1 ~ me ia ^ T3 )ip , (8) 

where a{z) parameterizes the kink, z is the spatial coordinate, and Tj denote 
the Pauli matrices. In this example, the induced fermion charge B^ of the 
ground state is given by the net change of a(z): 

BP) = / ^ dz = _ , (9) 

where Aa = a(+oo) — a(— 00). 

In a similar way, domain walls in a 3 + 1-dimensional theory can acquire 
a fermionic charge. Because of the planar symmetry, the computation can 
be reduced to that of the above 1 + 1-dimensional model. The starting point 
is the following simplified Lagrangian for the nucleon TV interacting with the 
non-fluctuating chiral field U: 

U = NidrfN ~ m N N L UN R - m N N R U+N L - \{N L N R ){N R N L ) , (10) 

where tojv is the nucleon mass. The last term is a four-fermion interaction 
term, and A > corresponds to repulsion in the U(l) channel. The 3 + 1- 
dimensional charge is given by 

s (4) = f N lQN d 3 x . (11) 



In order to reduce the computation of the four-dimensional charge 
to the two-dimensional problem of and (|^), we write the four-dimensional 
spinors N R and Nl in terms of a set of two-component spinors. Note that 
unless A 7^ 0, the contributions to i^iLirom different two-component spinors 
cancel (details will be given elsewhereOj). 

For simplicity, we consider a wall separating the true vacuum from its 
neighboring meta-stable ground state with —iTrlogU > 0. In the following 
section we will discuss under which circumstances these walls dominate over 
all other possible walls. 

Because of the coupling of ip to U, the effective nucleon mass ra e // takes 
on different values in different vacua k: 

m eff(k) = m N + m q f(6,k), (12) 
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where thc-function f(8, k) depends on the precise form of the effective potential 
V e ff(U)cB. In the adiabatic approximation, m e ff should be considered as a 
slowly varying function of z (the direction orthogonal to the domain wall). 
Fortunately, only the asymptotic values of m e ff enter the final answer. The 
resulting expression for the two-dimensional baryon number of the wall is 

= larccos^Hzt™, (13) 

where A has dimension of mass and can be found in terms of A and rriN- For 
the domain wall we are considering, is negative if we take the false vacuum 
to be at z = +00 and the true vacuum at z = —00. 

To find the original four-dimensional baryon charge B^ 4 \ we should take 
account of the degeneracy related to the symmetry under shifts along the wall 
plane. In many body physics the definition of the charge is B = J N l ^QN l d 3 x 
where the sum runs over all particles (possible quantum states). In our specific 
case this summation leads to the result 

flW = B^g { dxd V d P* d Pv = B W N , (U) 

J (27T)^ 

where g = 4 describes the degeneracy in spin and isospin. In the following, we 
shall estimate the value of N. 

For a wall with surface area S and for a fixed number of quantum states 
N we have 

AT Q f 9 S PF fTK\ 

^ = ^7(2^ = "ST' (15) 

where pp is the Fermi momentum. The Fermi energy Ep of the domain- wall 
fermions is determined by 

- f pd 2 p 2 4 [¥N 3 / 2 

The total energy of the fermions residing on the surface S is given by 

4 fW N 3 / 2 

E ° = ° s+ - 3 \h^w (17) 

The size of the surface which can accommodate the fixed number of fermions 
N can be found from the minimization equation 

—jJr\N=const — 0. (18) 
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which relates the density of fermions per unit area n = £ to the wall tension 
a: 

,3/2 



n 3/* = CT /rii. (19) 



Hence^the induced fractional charge on the domain wall follows from Eqs. (14) 
and @: 

Q = fiW = -|B( 2 )|7V = -SlB^la^A 1 / 3 , (20) 

47T 

which can be expressed in terms of a dimensionless constant ot\: 

Q = -SA% CD a u a, = -^-A^\Bi% (21) 



4 B-Shell Formation and Baryogenesis 

We now turn to our proposed baryogenesis (charge separation) scenario. At the 
QCD phase transition at the temperature T c ~ A-qcd the chiral condensate 
U forms. Because of the presence of nearly degenerate states, a network of 

Smain walls will arise immediately after T c by the usual Kibble mechanism 
At T c , the energy difference between vacua is negligible compared to the 
energy in thermal fluctuations, and hence the states are equidistributcd. The 
initial wall separation (correlation length £) depends on the details of the 
chiral phase transition but is expected to be microscopic. Below T Cl the wall 
network coarsens. We will assume that an infinite wall network will exist 
until a temperature at which time the energy difference between correlation 
volumes of the true vacuum and the false vacuum closest in energy becomes 
thermodynamically important. At this time, the wall network will decay into 
a number of finite clusters of the false vacuum which we will call B-shells. 

If 9 = 0, there are two degenerate meta-stable states \B > and \C > above 
the true vacuum of lowest energy \A >. For simplicity, we ignore meta-stable 
states of higher energy. A CP transformation exchanges the states \B > and 
\C >. Hence, the baryon charge of a A — B wall will be opposite of that of a 
A — C wall, and no baryon number will be left behind in the bulk because the 
number of A — B walls will be the same as the number of A — C walls. 

However, if at the temperature T c the value of 9 is different from 0, then 
the situation is very different. This is the case which will be considered below. 
Thus, we are making the assumption that the strong CP problem is cured by 
an axion at a temperature below T c . At T = T c , the axion is not yet in its 
ground state, and thus 9(T C ) might be of order unity. Note that as long as 
the initial value 9{T C ) is the same in the entire observed Universe, the sign of 
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the baryon asymmetry will also be the same. This will occur if the Universe 
undergoes inflation either after or during the Peccci-Quinn symmetry breaking. 
In this case there is a splitting of 8(T c )m q h^ CD between the energy densities 
of the states \B > and \C > which translates into a splitting AM ~ 6(T C )M 
between the masses of B-shells of the phases \B > and \C > with negative 
and positive baryon numbers (here, M stands for the B-shell mass at 8 = 0). 
We will assume that AM is larger than Td- Since the correlation length £ 
grows rapidly after T c , this requirement can be achieved without requiring a 
large value of 8. In this case, at the temperature T<j, only B-shells of one 
type, of negative baryon number, will remain. For a value £(T<j) ~ 10 T~ (a 
value which we argue below is reasonable) and assuming spherically symmetric 
B-shells, the criterium for 8(T C ) becomes: 



where in the last step we have replaced Td by T c to obtain a conservative 
bound. 

Note that the typical wall separation £(T<j) is rather uncertain since it 
depends on the initial correlation length at formation, on the details of the 
damping mechanism and on-^hCpiiitcrplay between the energy bias and the 
surface tension in the walls E3'E3'Eil. Given the typical wall separation £, the 
total area S in walls at the temperature Td within some reference volume V is 
S ~ y/f. The total baryon charge is given by ([H|). Since the entropy density 
is s = g*T^, where g* ~ 10 is the number of spin degrees of freedom in the 
radiation bath at Td, the net baryon to entropy ratio at Td becomes 



The evolution of the B-shells after Td requires a detailed study. Qualita- 
tively, we expect that the bubbles will shrink, but not decay completely since 
they will eventually be stabilized by the fermions. We expect the annihilation 
cross-section between a baryon and a B-shcll to be suppressed by a large power 
of the ratio of the Compton wavelength of the baryon and the radius of the 
B-shell. As the non-relativistic baryons can hardly cross the wall, we expect 
the shells to be stable against the escape of baryons from the interior, but able 
to lose heat by baryon pair annihilation and emission of the photons and/or 

Hutrinos. The quantum stability of the B-shells will be addressed separately 
, Generally, one expects an exponential suppression of quantum decays by 
the baryon charge of the surface. 




(22) 




(23) 
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To proceed, we introduce two dimensionless constants ai and «3, param- 
eterizing the total area and volume of the B-shells as a\V/t; and a^V, respec- 
tively. This parameterization does not imply any specific assumption about 
the form of the B-shells. Neglecting the expansion of the Universe between T c 
and T4 we obtain 

-i^~^ki- (24) 

Since the energy density ps in B-shells will red-shift as matter, the contribution 
£Ib of B-shells to the dark matter of the Universe is given by 

n PBjteq) PBjteq) rot -x 

Pr\teq) 9* d e 1 

where t eq is the time of equal matter and radiation. Assuming that the B-shell 
energy is dominated by the false vacuum energy, we obtain: 

S7 B ~a 3 ^f. (26) 

Comparing ( p4| ) and (p6|), we see that the resulting values of Ub/s and of VLb 
are related with each other via the geometric parameters «2 and a^. At the 
moment, we are not able to calculate these parameters directly. However, we 
can reverse the argument and ask what values of ct2 and are required in 
order to explain both £Ib ~ 1 and ub/s ~ 10 -10 . Taking Td ~ T c wc obtain 
03 ~ 10~ 6 and a 2 ~ 10 _6 ^T,j. Since by definition «2 and since £Td > 1, we 
are left with the window 

T- 1 < £ < lO 6 ^ 1 (27) 

for the proposed mechani sm tn be operative. This window is consistent with 
the Kibble-Zurek scenarioE3'E3 of defect formation in the early Universe. 



5 Discussion 



We have seen that it is possible, without fine tuning of parameters, to obtain a 
reasonable value of the baryon to entropy ratio in the bulk. B-shells will con- 
tribute to the dark matter of the Universe, and there is a region of parameter 
space for which B-shells will make up the bulk of the dark matter. Note that 
in our charge separation scenario, charges are separated only over microscopic 
scales. 

For the scenario to work, it is important that the B-shells be stable. We 
argue that Fermi pressure will stabilize the shells against collapse. We assume 
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that eventually the radius will be much larger than A.qq D , thus justifying the 
use of the results for the charge per unit area derived for a flat surface. In this 
case, the total energy for a spherical shell of radius R is 

(see ([l7])). Here, N is the total baryon number of the shell which is held fixed 
and can be estimated from the initial shell radius i?o = £, N = 47r£ 2 n, where 
the baryon number density n is given by (p^). By minimizing ( ffijj ) with respect 
to R (and noting that the surface tension term is negligible compared to the 
other two terms) , we can find the stabilization radius for a fixed initial baryon 
charge 

iV 3 / 2 

which then also determines the total energy of a B-shell, with the result that 
E ~ N g / S . Taking the value of £ to be at the upper end of the range ( p7|) we 
find E ~ 10 14 GeV or about 10- 10 g. 

Since the negative baryon charge is trapped in a topological configuration, 
we expect the annihilation cross section between nucleons and B-shclls to be 
greatly suppressed because of the mismatch between the Compton wavelength 
of the nucleons and the B-shell radius. Similarly, any charge loss mechanism 
will also experience this phase space suppression. Constraints on our proposed 
mechanism based on elastic scattering of B-shells in dark matter detectors 
remain to be explored. 

In conclusion, qualitative as our arguments are, they suggest that baryo- 
genesis can proceed at the QCD scale, and might be tightly connected with 
the origin of the dark matter in the Universe. 
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